NASA-CR-204774 


SwRI Proposal No. 15-7238 
February 27, 1997 




STUDIES OF DISKS AROUND THE SUN AND 

OTHER STARS 

NAGW-4468 


Submitted to: 

Dr. Joseph M. Boyce 
Code SR 

NASA HEADQUARTERS 
300 E. Street, S.W. 
Washington, DC 20546 


Prepared by: 

S. Alan Stern 

Southwest Research Institute 

Division 1 5 Boulder Extension Office 

Geophysical, Astrophysical, and Planetary Sciences 

1050 Walnut Street, Suite 426 

Boulder, Colorado 80302 

Telephone 1-303/546-9670 FAX 1-303/546-9687 



SOUTHWEST RESEARCH INSTITUTE 
Instrumentation and Space Research Division 
6220 Culebra Road, San Antonio, Texas 78238 
(210) 684-51 1 1 • FAX (210)647-4325 




INTRODUCTION 


This is the Final Report for NAGW-4468 (SwRI Project 15-7238), Studies of Disks Around 
the Sun and Other Stars, (S.A. Stern, PI). This is a NASA Origins of Solar Systems research 
program, and this NASA Headquarters grant has now been transferred to a new grant at NASA 
GSFC (NAG5-4082). Thus the need for this ” Final Report” on a project that is not, in fact, 
complete. 

We are conducting research designed to enhance our understanding of the evolution and 
detectability of comet clouds and disks. This area holds promise for also improving our under- 
standing of outer solar system formation, the bombardment history of the planets, the transport 
of volatiles and organics from the outer solar system to the inner planets, and to the ultimate fate 
of comet clouds around the Sun and other stars. According to “standard” theory, both the Kuiper 
Belt and the Oort Cloud are (at least in part) natural products of the planetary accumulation 
stage of solar system formation. One expects such assemblages to be a common attribute of other 
solar systems. 

Our program consists modeling collisions in the Kuiper Belt and the dust disks around 
other stars. The modeling effort focuses on moving from our simple, first-generation, Kuiper 
Belt collision rate model, to a time-dependent, second-generation model that incorporates physical 
collisions, velocity evolution, dynamical erosion, and various dust transport mechanisms. This 
second generation model is to be used to study the evolution of surface mass density and the 
object-size spectrum in the disk. 


PROGRESS 

1) We have now completed the first model of collision rates of the Kuiper Belt (KB). With this 
model we explored the rate of collisions among bodies in the present-day Kuiper Belt as a function 
of the total mass and population size structure of the Belt. We find that collisional evolution is an 
important evolutionary process in the KB as a whole, and indeed, that it is likely the dominant 
evolutionary process beyond ~ 42 AU, where dynamical instability timescales exceed the age of 
the solar system. Two further findings we report from this modelling work are: (i) That unless the 
Belt’s population structure is sharply truncated for radii smaller than ~l-2 km, collisions between 
comets and smaller debris are occurring so frequently in the Belt, and with high enough velocities, 
that the small body (i.e., km-class object) population in the disk has probably developed into 
a collisional cascade, thereby implying that the Kuiper Belt comets may not all be primordial, 
and (ii) that the rate of collisions of smaller bodies with larger 100 < R < 400 km objects (like 
1992QBi and its cohorts) is so low that there appears to be a dilemma in explaining how QB x s 
could have grown by binary accretion in the disk as we know it. Given these findings, it appears 
that either the present-day paradigm for the formation of the KB is failed in some fundamental 
respect, or that the present-day disk is no longer representative of the ancient structure from 
which it evolved. In particular, it appears that the 30-50 AU region of the Kuiper Belt has very 


1 




likely experienced a strong decrease in its surface mass density over time. This in turn suggests 
the intriguing possibility that the present-day Kuiper Belt evolved through a more erosional stage 
reminiscent of the disks around the A-stars (3 Pictorus, a PsA, and a Lyr. These results were 
published in The Astronomical Journal in 1995 and 1996. 

2) We have used our intial collision rate model and a second code to estimate the detectability of 
IR emission from debris created by collisions. We found that eccentricities in the Kuiper Belt are 
high enough to promote erosion on virtually all objects up to ~ 30 km, independent of their impact 
strength. Larger objects, such as the 50-170 km radius “QBi” population, will suffer net erosion if 
their orbital eccentricity is greater than « 0.05 (~ 0.1) if they are structurally weak (strong). The 
model predicts a net collisional erosion rate from all objects out to 50 AU ranging from 3 x 10 1 
to 10 19 g yr, -1 depending on the mass, population structure, and mechanical properties of the 
objects in the Belt. We find two kinds of collisional signatures that this debris should generate. 
First, there should be a relatively smooth, quasi-steady-state, longitudinally isotropic, far IR (i.e, 
~ 60 fim peak) emission near the ecliptic in the solar system’s invariable plane ecliptic, caused by 
debris created by the ensemble of ancient collisions. The predicted optical depth of this emission 
could be as low as 7 x 10~ 8 , but is most likely between 3 x 10 7 and 5 x 10 3 * * 6 . We find that this 
signature was most likely below IRAS detection limits, but that it should be detectable by both ISO 
and SIRTF. Secondly, recent impacts in the KB should produce short-lived, discrete clouds with 
significantly enhanced, localized IR emission signatures superimposed on the smooth, invariable 
plane emission. These discrete clouds should have angular diameters up to 0.2 deg, and annual 
parallaxes up to 2.6 deg. Individual expanding clouds (or trails) should show significant temporal 
evolution over timescales of a few years. As few as zero or as many as several 10 2 such clouds 
may be detectable in a complete ecliptic survey at ISO’s sensitivity, depending on the population 
structure of the Kuiper Belt. This work was published in Astronomy & Astrophysics in 1996. 

3) We then employed our model to study the collisional environment in the ancient Kuiper Belt. 

We explored the consequences of a massive, primordial Kuiper Belt using a collision rate model 

that assumes the dominant growth mechanism in the 35-50 AU region was pairwise accretion. We 

found that the growth of QBi -class objects from seeds only kilometers in diameter required a very 

low eccentricity environment, with mean random eccentricities of order 1% or less. Duncan et al. 
(1995) have shown that the presence of Neptune induces characteristic eccentricities throughout 

the 30-50 AU region of a few percent or greater. We therefore concluded that growth of objects in 
the 30 to 50 AU zone to a least this size must have occurred before Neptune reached a fraction of its 
final mass. Once Neptune grew sufficiently to induce eccentricities exceeding ~1%, we found that 
the disk environment became highly erosive for objects with radii smaller than ~20— 30 kilometers, 
which likely created a flattening in the disk’s population power law slope between radius scales of 
«30 to «100 km, depending on the density and strength of such objects. This erosive environment 
could have resulted in sufficient mass depletion to evolve the disk to its present, low- mass state, 
independent of dynamical losses (which surely also played an important role). During the period 
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of rapid erosive mass loss, the disk probably exhibited optical depths of 10 4 to 10 5 (reminiscent 
of P Pictoris), for a timescale of ~10 7 to ~10 8 years. As a result of the evolution of the disk inside 
50 AU, we suggested that (i) the present-day Solar System’s surface mass density edge near 30 
AU is actually only the inner edge of a surface mass density trough., and (ii) that the surface mass 
density of solids may rise back beyond ~50 AU, where the giant planets have never induced erosive 
high eccentricities. Indeed, the growth of objects in the region beyond 50 AU may be continuing 
to the present. This work was published in The Astronomical Journal in late 1996. 

4) In collaboration with Col Colwell, I have completed the main task of the proposed work in this 
grant effort, to construct a time-dependent code for modelling of Kuiper Belt collisional evolution. 
Two papers are now in press discussing key results obtained with this model (Stern k. Colwell 
1997a, b). Applying a time-dependent model of collisional evolution of the EKB, we found that 
under a wide range of assumptions, collisional evolution should have depleted the mass of the 
30-50 AU zone by >90% early in the history of the solar system, thereby creating a deep scar or 
gap in the surface mass density across a wide region beyond Neptune, much like what is observed 
today. Dynamical erosion may have further accelerated the depletion process. Given the fact 
that Neptune has had far less dynamical influence beyond 50 AU, our results also suggest that 
unless the solar nebula was truncated near 50 AU, then surface mass density of solids beyond ~50 
AU increases again, most likely dramatically. In paper II, we employed the new, time-dependent 
collisional evolution code to study the conditions under which the ~50— 200 km radius Edgeworth- 
Kuiper Objects (EKOs) in the region between 30 and 50 AU (now called the Edgeworth-Kuiper 
Belt, or EKB; Edgeworth 1943, 1949; Kuiper 1951) were formed. Assuming that these bodies were 
created by pairwise accretion from 1 to 10 km building blocks, we find that three conditions were 
required, namely: (i) at least 10 and more likely 35 of solids in the primordial 30 to oO 
AU zone, (ii) mean random orbital eccentricities of order 0.002 or smaller, and (iii) mechanically 
strong building blocks. Furthermore, we find that the accretion of 100-200 km radius bodies in 
the 30 to 50 AU region from collisions among a starting population of 1 to 10 km building blocks 
required «10 8 -10 9 years, with the lower range only being reached for 30 to 50 AU zone masses 
approaching 100 of solids or mean random orbital eccentricities <0.005 (which we do not 
believe is realistic after gas dissipation). Therefore, unless accretion had already produced many 
building blocks significantly larger than 10 km in diameter at the time the nebular gas was removed, 
our results also indicate that Neptune did not form on a timescale much shorter than ~70 Myr, 
and could well have required many hundreds of Myr to approach its final mass. We also explore 
the growth of Pluto-scale (i.e. , radius 1000-1200 km) objects in the 30 to 50 AU region under a 
variety of assumptions. We further find that once ^300 hundred kilometer radius objects were 
formed, the growth of 1000 km radius and larger objects occurs relatively easily and comparatively 
quickly. The lack of many Plutos in the 30 to 50 AU zone therefore argues strongly that growth was 
terminated in that region rather abruptly at the time the presently observed population of 100-200 
km radius EKOs were being completed. In the region beyond 50 AU where Neptune’s dynamical 
influence was much reduced, model runs yield 100 km to 1000 km radius, and perhaps even larger 
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bodies which should be detectable with on-going or soon-to-be started surveys. We suggest that if 
dynamical conditions did not remain calm enough to allow Pluto itself to be grown in the 30 to 50 
AU zone before perturbations from Neptune created a dynamically erosive, low-mass environment 
there, then it may be that Pluto was grown beyond the influence of Neptune s perturbations and 
later transported inward, perhaps in part via the Charon-forming collision. 

5) These five papers were accompanied by an invited review, submitted to the Planetary Ices book, 
summarizing the present state of knowledge about the Kuiper Belt and Pluto, and another review, 
on the origin of Pluto, which is now in press for the UA Space Science serives volume, Pluto & 
Charon .. 

6) We organized and sponsored a 2-day workshop on collisions in the Kuiper Belt. This workshop 
was attended by D. Davis (PSI), P. Farinella (Italy), R. Canup (U. Colorado), M. Festou (France), 
J. Colwell (U. Colorado), H. Levison (SwRI), and PI Stern (SwRI). The proceedings of this work- 
shop were informally published and the distributed among the participants. A copy was also sent 
to Origins program scientist Trish Rogers. 

7) We have also written a popular-level article on extra-solar comets for Astronomy magazine. 

8) Additionally, PI Stern has given nine invited talks summarizing the collisional modelling results 
obtained under the Origins program. A list of these invited talks is attached. 

9) For the remainder of this program we plan to (i) concentrate on improving our KB collisional 
evolution code to include coupled mass-velocity evolution, and (ii) to then begin exploiting the 
improved code to better understand the growth of objects in the Kuiper Belt. 
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ABSTRACT 

We explore the rate of collisions among bodies in the present-day Kuiper Disk as a function the total mass 
and population size structure of the disk. We find that collisional evolution is an important evolutionary 
process in the disk as a whole, and indeed, that it is likely the dominant evolutionary process beyond *-42 
AU, where dynamical instability time scales exceed the age of the solar system. Two key findings we report 
from this modeling work are: (i) That unless the disk’s population structure is sharply truncated for radii 
smaller than ~ 1-2 km, collisions between comets and smaller debris are occurring so frequently in the disk, 
and with high enough velocities, that the small body (i.e., KM-class object) population in the disk has 
probably developed into a collisional cascade, thereby implying that the Kuiper Disk comets may not all be 
primordial, and (ii) that the rate of collisions of smaller bodies with larger 100</?<400 km objects (like 
1992QB! and its cohorts) is so low that there appears to be a dilemma in explaining how QBj’s could have 
grown by binary accretion in the disk as we know it. Given these findings, it appears that either the 
present-day paradigm for the formation of Kuiper Disk is failed in some fundamental respect, or that the 
present-day disk is no longer representative of the ancient structure from which it evolved. This in turn 
suggests the intriguing possibility that the present-day Kuiper Disk evolved through a more erosional stage 
reminiscent of the disks around the stars /? Pictorus, a PsA, and a Lyr. © 1995 American Astronomical 
Society. 


1. INTRODUCTION 

Over the past few years, both the theoretical underpin- 
nings and the observational evidence for a disk of comets 
and larger bodies beyond the orbit of Neptune has become 
increasingly secure (Jewitt & Luu 1995; Cochran et al 
1995). It now appears assured that the solar system possesses 
such a disk of material, and that this region is likely to con- 
tain the source population for the low-inclination, short- 
period, Jupiter Family Comets (Stem 1995a). 

In this paper I explore the rate at which objects collide in 
the Kuiper Disk region. The basic rationale for such a study 
is rooted in the combination of a K^-IO 4 times higher num- 
ber density of comets and 10 1 times average orbital speed in 
the KD, compared to the Oort Cloud (Stem 1988), which 
together imply that collision rates should be 10 7 -10 9 times 
higher in the Kuiper Disk. Further rationale is provided by 
analogy to the asteroid belt. The average surface mass den- 
sity in the Kuiper Disk ( — 10 23 — 10 24 g AU" 2 ) is similar to 
the value of —BXIO 23 gAU” 2 in the asteroid belt, where 
collisions play an important and well-known evolutionary 
role. Even accounting for the ~4 times lower random veloci- 
ties at 40 AU in the Kuiper Disk (as opposed to 2 AU in the 
asteroid belt), the collisional intensity in the Disk (i.e., col- 
lisions cm” 2 s” 1 on a given target) is not very different from 
the asteroid belt. 

Among the questions about the Kuiper Disk that one 
wishes to address with collision rate modeling are: What is 
the rate of collisions in the disk today? Is the Kuiper Disk 
collisionally evolved; that is, are cratering collisions an im- 


portant surface modification process in the disk, and is the 
rate of collisions high enough to permit evolution in the size 
spectrum of bodies in the disk? Is it possible to constrain the 
properties of the ancient disk via collisional results? Is it 
possible to constrain the properties of the distant, as-yet un- 
detected reaches of the disk via collisional results? And, are 
there detectable signatures of these collisions? 

This paper represents an initial attack on several of these 
questions. It is organized as follows: In Sec. 2, I briefly re- 
view the evidence for the Kuiper Disk; Section 3 describes a 
model for computing collision rates in the present-day 
Kuiper Disk; Section 4 describes the results of model runs 
for the present-day Kuiper Disk; Section 5 examines the im- 
plications of these results; Section 6 explores whether colli- 
sions in the present-day disk promote accretion or erosion; 
Section 7 summarizes the results obtained in this paper and 
points out two significant inconsistencies between the colli- 
sional modeling results obtained here and the present under- 
standing of the origin of objects in Kuiper Disk. Among the 
implications of the work reported here is that the present-day 
disk appears to be the remnant of a former disk with more 
mass, and very likely lower mean eccentricities, than ob- 
served today. 

2. THE KUIPER DISK 

Almost a half-century ago, Edgeworth (1949), and later 
Kuiper (1951), made prescient predictions that the Sun 
should be surrounded by a disklike ensemble of comets and 
other “debris” located beyond the orbit of Neptune. The case 
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for such a primordial reservoir was strengthened when it was 
pointed out that such a disk could be an efficient source 
region to populate the low-inclination, short-period comets 
(Fernandez 1980). Convincing dynamical simulations sup- 
porting this link between most short-period comets, and the 
Kuiper Disk (KD) region, first appeared when Duncan et al. 
(1988) and later Quinn et aL (1991) showed that a low- 
inclination source region appears to be required for the low- 
inclination orbit distribution of the Jupiter Family Comets. 
Figure 1 is a schematic diagram revealing the gross architec- 
ture of the disk in relationship to the orbits of the five known 
outer planets. 

The computational capabilities available to Duncan and 
co-workers in the mid-1980s required some important ap- 
proximations be accepted (for reviews of this work, cf. 
Weissman 1993, and Stem 1995a). These compromises were 
criticized by Bailey & Stagg (1990), but Duncan et al .' s 
work generated interest in the Kuiper Disk by both modelers 
and observers. Of particular relevance are the Holman & 
Wisdom (1993) and Levison & Duncan (1993) studies of 
orbital evolution in the disk. These two groups found a time- 
dependent dynamical erosion of the disk population inside 
~42 AU, caused by nonlinear perturbations from the giant 
planets. The dynamical chaos resulting from these perturba- 
tions is ultimately responsible for the transport of short- 
period comets from the long-lived Kuiper Disk reservoir to 
planet crossing orbits where they can be routinely detected. 
Based on the bias-corrected population of Jupiter Family 
comets and the dynamical transport efficiency of comets 
from the Kuiper Disk to the inner planets region, Duncan 
et aL (1995) have estimated that 6X 10 9 comets orbit in the 
Disk between 30 and 50 AU from the Sun. 

Observational confirmation of the Kuiper Disk was first 
achieved with the discovery of object 1992QBj by Jewitt & 
Luu (1993). As of early 1994, no fewer than 25 QB r like, 
frans-Neptunian objects have been discovered (Jewitt & Luu 
1995; Stem 1995a). These icy outer solar system bodies are 
expected to have dark surfaces consisting of an icy matrix 
contaminated by silicates and organics. Assuming a typical 
(i.e., cometary) geometric albedo of 4%, and the absence of 
coma, the distances and magnitudes of these objects indicate 
they have radii between roughly 50 and 180 km. Based on 
the detection statistics obtained to date, one can easily esti- 
mate that a complete ecliptic survey would reveal — 3.5 X 10 4 
such bodies orbiting between ^30 and 50 AU. Simple 
power-law extensions of this population predict a cometary 
population (which we define as bodies with radu between 1 
and 6 km) of — 10 10 , which is similar to dynamical modeling 
results obtained by Duncan et aL (1995) to satisfy the short- 
period comet flux. Very recently, Cochran et aL (1995) have 
reported Hubble Space Telescope results giving the first di- 
rect evidence for comets in the Kuiper Disk. 


3. COLLISION RATE MODELING 

Our model for estimating collision rates in the Kuiper 
Disk begins by defining the Disk in terms of a power law 
exponent, a, on the size distribution of objects in the disk, so 



Fig. 1. Schematic depiction of the Kuiper Disk and the orbits of the outer 
planets, including Pluto. The clearing between the orbit of Neptune and the 
inner edge of the present day disk is created by the dynamical perturbations 
of the giant planets (cf. Holman & Wisdom 1993; Levison & Duncan 1993). 
The position of the outer boundary of the Kuiper Disk is not well con- 
strained, and may well extend much farther than shown here. 

that the number of objects dN(r ) between radius r and 
r+dr is given by 

dN(r)=N 0 r a dr, 0) 

where V 0 is a normalization constant set by the estimated 
number of QB! objects. For the runs presented here, the bod- 
ies in each successive size bin r are a factor of 1.6 times 
larger in size (and equivalently, 4 times higher in mass). We 
a priori assume a size range beginning at r=0.1 km, and 
extending upward to r= 162 km. 

We also define a power law exponent {3 on the radial 
distribution of surface mass density 2(/?) in the disk, so that 

2(/0=Xo**, (2) 

where 2o is the normalization constant. 

Once an input disk is defined as described above, the 
model bins the disk into a series of concentric tori that are 1 
AU in radial width. For each radius bin/heliocentric bin pair, 
the model computes the collision rate a target will experience 


Table 1. Collision model run cases. 


Model 

30 AU < R <50 AU 

Tw 

Di* 

Ttp« 

Nq»i« 

30 AU < R <60 AU 

* 

30 AU < A <50 AU 

CSD let IS 

0.16 U* 

NOM 

CUB 

30,401 

txl* 

CSD let 26 

0.12 U* 

NOM 

DMB 

27,740 

TxlOP 

CSD let 13 

0.42 bU 

CM 

CUB 

41,102 

ix 10“ 

CSD let 23 

042 M* 

CM 

DMB 

31416 

3 x 10** 

CSD 3e4 15 

0 07 U* 

NOM 

CUB 

17,950 

4x 10* 

CSD lefl 13 



CM 

CMB 

1.2 x 10* 

14 x 10 11 

Notes to Table 1 


^S Aiak is the integral mass over all size bins. CMB= constant mass per radial 
bin (/?= - 1); DMB= constant mass per radial bin (/?= -2). NOM= Nominal 
mass per size bin (<*=-ll/3); CM = constant mass per size bin (ar=-4). 
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are for a disk with (e)**0.20. 


from objects in all bins of equal or smaller size. This is of 
course a function of (e), since (e) controls both the internal 
velocity dispersion in the disk, as well as the degree of he- 
liocentric bin crossing. In what follows we assume (i) 
as %(e). Because there is presently no information on the way 
in which ensemble-averaged inclinations ({i)) and eccen- 
tricities ((e)) vary in the Kuiper Disk, we adopt a disk-wide 
(i) and (e) for each run, and vary these quantities from run 
to run as free parameters to explore how sensitive the model 
results are to these variables. 

Collisions are not allowed outside the boundaries of the 
disk, so in the case of moderate or high eccentricity orbits, 
objects can spend significant time in “open” space outside 
the disk where collisions are not allowed to occur. This cre- 
ates edge effects, but such effects may actually occur if the 


disk in fact sharply truncates at its boundaries. 

To compute collision rates we adopt a particle in a box 
formalism. In this approach, the instantaneous collision rate 
c of objects with semimajor axis a, eccentricity e , and ra- 
dius r x being struck by objects of radius r y is 

c(r x ,r y ,a,e) = ncr g v, (3a) 

where n is the local space density of impactors, v, is the 
local average crossing velocity of the target body against the 
KD population at distance R, and cr g is the collision cross 
section of the impactor+ target pair, corrected for gravita- 
tional focusing. Gravitational focusing is an important cor- 
rection for targets in the QB, size range and larger, particu- 
larly in the case of very low disk eccentricities (e.g., 
(e)<10~ 2 ). The orbit-averaged collision rate c(r x ,r y ,a,e) 
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can be written to show its implicit dependencies in the model 
as: 


«<!+<*» 

c(r x ,r y ,a y e) = 2 /(a,(e),/?)n(r > ,/?) 

xy »^CSC(X + y))’ 

(3b) 

Here the term / represents the ratio of T(a f (e),R)> the time 
the target body in an orbit defined by (a,{e),{/)) spends in 
each torus it crosses during its orbit, to this target’s orbital 
period, (47r 2 a 3 /GAf 0 ) 1/2 . 

To compute T(a y (e) y R y I solve the central-field Kepler 
time of flight equation explicitly for every (a y (e)) pair in the 
run parameter space. The number density of impactors 
n(r y ,/?) in the torus centered at distance R is computed from 
the defining mass of the disk, its wedge angle (i), its radial 
surface density power law, and the population size structure 
power law. 

To compute the average crossing velocity of the impactor 
population on the target body when the target is in the bin at 
heliocentric distance /?, we use 

v xy =v K (a ) > / 2 <*> 2 + 2 </> 2 - 3 ( ^f] 7 (4) 


(Petit & Henon 1987), where v K is the average Keplerian 
orbital speed of the target body, and the term under the radi- 
cal is the relative velocity correction for crossing orbits. The 
collision cross section cr g is computed according to 


<r t =Tr(r x +r y ) 2 


1 + 


2G(/n x + m y ) 
v^(r^+r~) 


( 5 ) 


where the term in brackets adjusts for gravitational focusing. 
To compute masses from radii I assume a density of 1 
g cm -3 . 

As a result of these calculations and the nested loops in 
a y r x , and r y , the model produces an array of collision rates 
c(r x ,r y ,a,e) throughout the specified disk, where the free 
parameters defining the disk are the total number of QB^s 
interior to 50 AU, a, 0, and (e). From this array the model 
computes subsidiary quantities such as the mean time be- 
tween collisions t < r x ,r y ,a,e), 

r(r x ,r y ,a,e) c (r x ,/‘ y ,n,e), ^ 

the mass impact rate from all impactors on each target size 
class: 


M(r x ,a,e)= 2 c(r x ,r y ,a,e)m(r y ), (7) 

and the total collision rate on entire population in each target 
size bin: 

C(r x ,a,e)= 2 c(r x ,r y ,a,e)N{r x ,a), (8) 

r »* r « 

where N(r x ,a) is the population of targets of radius r x with 
semimajor axis a. We also compute a characteristic time for 
growth, t g , as 


T G (r x ,d,e) = 


Af(r x ,a) 

7}M(r x ,a,e)' 


where tj is the mass accretion efficiency per collision. 


(9) 


4. MODEL INPUT PARAMETERS 

In the runs presented below, I assume a disk inner radius 
of 35 AU and an outer radius of 70 AU. I let eccentricity 
range as a free parameter from 1X10 4 to 2X10 \ which 
extends over the range of detected eccentricites of QI^ ob- 
jects detected to date (Jewitt & Luu 1995; H. Levison, per- 
sonal communication 1995). As noted above, I assume the 
equilibrium condition ( f) = K e ) ■ 

Four cases defining the radial mass dependence and size 
distribution of objects in the disk have been studied. These 
four cases represent the various combinations of two radial 
mass distributions [cf., a in Eq. (1)] and two size distribu- 
tions [cf„ /3 in Eq. (2)] . 

For the radial distribution of mass in the disk, the two 
cases we run are defined as follows One case assumes a 
constant mass per radial bin (CMB; a— — 1), which corre- 
sponds to a surface mass density that declines with heliocen- 
tric distance as The second, and more realistic case, 
assumes a declining mass per bin (DMB; a =— 2), corre- 
sponding to a surface mass density falling like R 2 . These 
two cases bracket the realistic range of parameter space (Lis- 
sauer 1987). 

Concerning the size distribution of objects in the disk 
population, the model grid allows for 17 size bins, each a 
factor of 1.6 larger in radius. We assume a minimum radius 
for KD imp actors of 0.1 km. This results in an upper size 
limit of r=162 km, which is consistent with the largest de- 
tected bodies among the QB t population. Our favored size 
distribution, which we call the nominal (NOM) case, con- 
nects the observationally estimated 3.5 x 10 4 QB (-sized ob- 
jects (Jewitt & Luu 1995) inside 50 AU with the modeling- 
derived estimates of ~10 10 comets (Duncan et al. 1995) in a 
single power law with a=-ll/3. Our second case assumes 
a=- 4, which gives a constant mass in every logarithmic 
size bin; this case is called the CM case. Relative to the 
NOM case which produces **" 10 10 for 35,000 QB(’s (100 km 
in radius or larger), the CM case produces “*5X 10 10 comets. 

Table I summarizes some the important attributes of these 
four run cases, as well as two additional run cases described 
in Secs. 6 and 7. With these preliminaries described, we now 
discuss the results relating to these four model cases. 

5. MODEL RESULTS: COLLISION AND GROWTH TIME SCALES 

Figures 2 and 3 depict the collision time scale results 
obtained using the model described in Secs. 3 and 4. Results 
are presented at two heliocentric distances, 40 AU (on the 
left) and 60 AU (on the right). In each figure, the upper 
panels show the collision time scale for («)’'* 10 2 , and the 
lower panels show the collision time scale for {«)<** 2 X 10“ . 
These values of (e) bound the measured eccentricity of all 
Kuiper Disk objects with known eccentricities. Similar data 
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Fig. 3. Contours of the collision time < 
model run with a constant mass per heliocentric bin (CMB) and constant i 
the lower panels are for a disk with {e)~0.20. 


scales (in years) at two locations in the Kuiper Disk (40 and 60 AU). as a function of target and impactor size for the 
rf!nr«,«ie bin fCMB) and constant mass (CM) per bin size structure. The upper panels are for a disk with (e)-0.01; 


that have been computed for the two DMB 2 ) cases 

are not shown because the results are not significantly differ- 
ent 

These data can be used to ascertain a number of interest- 
ing facts about collisions in the present-day disk. Two results 
that are relevant to our later discussions concern the follow- 
ing. 

(1) Collisional Time scales on Comets in the Disk: We 
define a “comet” as those disk objects in the radii bins from 
1 to 6 km. In the case of the NOM population structure (cf., 
Fig. 2), the largest impactor a comet at 40 AU typically col- 
lides with in 4x 10^ yr has a radius times smaller than the 
comet itself; at 60 AU the largest impactor on a comet is 
typically — 10 times smaller. In the case of a CM population 
structure (cf., Fig. 3), which has more small bodies and 


therefore shorter collisional time scales than the nominal- 
case population, comets are expected to be struck by ap- 
proximately like-sized impactors at 40 AU, and 2.5 times 
smaller objects at 60 AU. And, 

(2) Collisional Time scales on QB l Bodies: We define 
“QBi bodies” to be objects in the 102 and 162 km radius 
bins, which span essentially the range of detected QB t radii 
(see, e.g., Jewitt & Luu 1995). Notice in Figs. 2 and 3 that 
over the age of the solar system, the largest impactor on a 
typical QBj body will be -6-16 km in radius, depending on 
the population structure and eccentricity of the disk. Notice 
also that each QBj object will suffer a cratering collision 
with a km-class object every 10 6 -10 7 yr in a CM population 
and every — 10 7 -10 8 yr in the NOM population. Among the 
entire population of ~3.5X10 4 QBj bodies inside ~50 AU, 
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and a nominal size structure. 


one expects ~ 10 2 — 10 3 collisions with a kin-class objects ev- 
ery year, depending on whether the population structure is 
more like the NOM or CM case. These collision rates sug- 
gest that although impacts on individual objects occur infre- 
quently, the population ensemble produces collisions fre- 
quently. This in turn suggests that a significant amount of 
dust may be injected into the Kuiper Disk every year, possi- 
bly leading to detectable signatures. This subject is beyond 


the scope of this paper, but is thoroughly investigated else- 
where (cf.. Stem 1995b). 

The results presented in Figs. 2 and 3 demonstrate very 
clearly that both small and large objects in the Kuiper Disk 
suffer collisions on time scales much shorter than the age of 
the solar system. 

It is next crucial to ask whether present-day rate of colli- 
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FIG. 5. Lower limit growth timescales as a function of target size as heliocentric distance for the model run with a Constant Mass per heliocentric bin (CMB) 
and constant mass (CM) per bin size structure. 


sions is large enough to have built the largest (i.e., QBi) 
bodies we see in the disk. To address this question I have 
computed growth times using the formalism imbedded in Eq. 
(9), with the assumption that the growth efficiency factor 
(i.e. T the mass accreted divided by mass incident) is unity. 
With 77 = 1 , collisions are completely inelastic. This is physi- 
cally unrealistic, since many collisions will result in erosion 
of the target rather than net accretion; however, it provides a 
useful lower limit to the actual growth times. As we shall see, 


even the lower limit QBt growth times are longer than the 
age of the solar system. 

Figures 4 and 5 depict the results of such lower-limit 
growth time calculations for the same two model runs that 
produced the collision time scales in Figs. 2 and 3, respec- 
tively. 

The results shown in Figs. 4 and 5 can be summarized as 
follows. For both the NOM and CM population size struc- 
tures, collisions are so infrequent that even km-scale bodies 






cannot not accrete their own mass in the age of the solar 
system, even if every collision is perfectly accretional (i.e., 
inelastic). In the case of the CM population structure, the 
largest objects that can be grown in 4X 10 9 yr are only a few 
km in radius. 

These results are not a strong function of (e) if { e ) 
>0.01, as appears to be the case in the present-day disk. As 
a result , we must conclude that binary accretion in the 
present day disk cannot explain the growth of QB r class bod- 
ies on time scales less than about an order of magnitude 
longer than the age of the solar system. The implications of 
this finding will be discussed in more detail in Sec. 7. 


Table 2. Critical eccentricities (e*) for erosion. 


Target Radius 

35 AU 
(Strong) 

80 AU 
{Strong^ 

35 AU 
(Weak) 

60 AU 
(We*k) 

001 km 

7 x 10-* 

9 x 10-* 

1 x 10-* 

1 x 10-* 

010 km 

6 X 10"* 

7 x 10“* 

2 x 10~* 

3 x 10-* 

100 km 

5 x 10-’ 

6 x 10-* 

2 x 10" 1 

3x10-’ 

170 km 

9 x 10- J 

Ijx 10-* 

4 x 10- J 

5 x 10-’ 


Notes to Table 2 

Strong implies p=2 gem' 3 and r=3xi0 6 erg g _1 ; weak implies p=0.5 
gem -3 and j= 3X10 4 erg g' 1 . In both cases we take *=8 and v^-O.lOv^ 
(e.g., Davis et ai 1989); see Sec. 6 for additional details. 


6. ECCENTRICITIES FOR GROWTH AND EROSION 

Up to this point we have not been strongly concerned with 
the issue of whether collisions in the KD promote commu- 
nition or growth. Instead we have been satisfied to simply 
count collisions and compute time scales. We have seen that 
binary collisions are too infrequent to explain the growth of 
objects larger than a few km in radius, even if all collisions 
promote growth. Now we explore whether growth can take 
place at all, or whether instead the collisions promote ero- 
sion. 

Whether a given collision between an impactor and a tar- 
get results in growth or erosion depends primarily on the 
energy of the impact In the Kuiper Disk a typical approach 
velocity of two objects at a distance large compared to the 
Hill sphere of the target can be reasonably -well approxi- 
mated by 

U .= V2M<* 2 )+<i 2 » ,/2 . (10) 

where v K is the local Keplerian velocity. For the standard 
assumption that (/) = K«)» we have, 

u.= V3<e 2 > , 'V (11) 

The energy at impact is therefore given by 

WL P =IAt(««*+^ 2 > W2 w*) 2 . (12) 

where p is the reduced mass and is the escape velocity 
of the two colliding bodies measured at the radius of impact. 
The critical velocity for net erosion to occur is given by the 
requirement that the specific impact energy must exceed the 
combined energy lost (a) to dissipation, (b) to break up the 
surface, and then (c) to disperse the ejecta out of the gravi- 
tational well of the combined mass of the impactor/target 
collision pair. The impact energy, as given by Eq. (12), 
must equal or exceed these energy sinks; if it does not, the 
target will accrete some mass in the collision. The critical 
condition for the target to lose mass occurs when the mass of 
the ejecta exceeds the mass of the impactor. Therefore, if the 
impactor mass is small compared to the target, we require 

+ y u^) f (13) 

where v s represents the velocity required to mechanically 
shatter the target surface, represents the velocity required 
to disperse the debris to infinity, and k is a factor that takes 
into account energy losses partitioned into heat, sublimation, 
hydrodynamic effects, and other factors. We take the specific 


energy for mechanical breakup of the target surface to be 

v s = 77 , (14) 


where s is the specific strength of the target material at zero 
compression. And of course the escape velocity is given by 


1 2GA/ t 
V r t 


(15) 


where G is the universal gravitational constant, M, is the 
combined mass of the target and impactor, and r, is the com- 
bined radii of these two objects. From Eqs. (9)— (15) one can 
derive the condition which we must solve for 

vl K +2yJ{3)e*VicV ac +3(e*) 2 v\-KV^-KV^-0, (16) 

to obtain e*, the critical eccentricity at which impact ener- 
gies are high enough to promote net erosion. Notice e* is a 
function of several parameters, including the target strength, 
size, and mass, as well as the heliocentric distance. 

Table 2 gives solutions to Eq. (16) for the critical erosion 
eccentricity e*, both for impacts onto strong (e.g., rock/ice) 
targets (p=2 gem -3 and r=3xi0 6 erg g -1 ), and relatively 
weak (e.g., snowlike) targets (p= 0.5 gem -3 and s=3X10 4 
erg g -1 ) at heliocentric distances of 35 and 60 All. Follow- 
ing the results discussed in Fujiwara et aL (1989), we assume 
u ej =0.2r; iinp and *=8. 

The results presented in Table 2 show that e* for 
QB[ -sized targets with radii near 100 km, £*2*0.02-0.03 is 
required for net erosion if they are weak, and 1 *^ 0 . 05 - 0.06 
is required if they are strong. Similarly, for QB,-like objects 
with /?= 170 km, which is comparable to the largest- 
discovered objects in the disk to date, £*2*0.04-0.05 is re- 
quired to result in net erosion if the objects are weak, and 
£*2*0.09-0.10 is required if the objects are strong. For ref- 
erence, at 35 AU an (£)=0.01 corresponds to a typical en- 
counter velocity at infinity of 87 m/s. We conclude from 
these results and the orbits of objects detected to date that 
some QBj’s should be undergoing erosion, while others may 
be in an accretional regime, depending on their eccentricity 
and strength. However, it is worthwhile to note that if the 
characteristic ejecta velocity v ej is as low as a few percent of 
the impact speeds v^, then e* will rise dramatically and the 
QB| population will be in an accretional mode, even for 
eccentricities as high as 0.2-0.5. Unfortunately, until much 
better eccentricity statistics become available, it is not pos- 
sible to determine if the QB| population as a whole is gain- 
ing or losing mass. All we can say is that the range of de- 
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tected eccentricities spans the range of e*’s, creating a 
complex situation. 

The results presented in Table 2 show that e* for comets 
is in the neighborhood of 10“ 3 to 1(T 2 , depending in large 
part on the true strength of comets. These results remain 
valid even if the characteristic ejecta velocity is as low as 5% 
of the impact speed, instead of 20%, as assumed in Table 2. 
This result and the fact that cometary orbit inclinations in the 
disk appear to be like QB! inclinations, imply that collisions 
on comets today arc erosive. This finding also indicates that 
a collisional cascade is probably taking place among the 
small bodies in the Kuiper Disk. As pointed out by P. 
Farinella (personal communication 1995), this finding 
strengthens the analogy made in Sec. 1 between the Kuiper 
Disk and asteroid belt collisional regimes. 

To determine how much mass a typical comet will loose 
in the age of the solar system, we combine the collision time 
scales in Figs. 2 and 3 with the algorithm outlined in Eqs. 
(10)— (16) to calculate a characteristic time scale (M/M) for 
such objects to erode to zero mass. This is accomplished 
through a numerical code, which we point out, only allows 
mass loss when e>e*. With this code, we find that between 
35 and 55 AU, the critical size for catastrophic (i.e., com- 
plete) erosion is —1-2 km, depending on the properties of 
the target and the disk population structure. In addition, we 
find that comets perhaps as large as 4 km in radius can ex- 
hibit erosion timescales shorter than the age of the solar sys- 
tem inside —40 AU, if (e)>0.04. 

To support these conclusions. Fig. 6 shows a set of ero- 
sion time scales calculated for a mechanically strong (i.e., 
j=3X10 6 erg g” 1 ) comet 1 km in radius, assuming 
u e j=0.1u iinp . The impact time scales used in this calculation 
were from the Fig. 2 dataset. Figure 6 shows that throughout 
the region from 35 to 55 AU, the erosion time for such bod- 
ies is less than or equal to the age of the solar system. These 
erosion time scales will be further shortened if either comets 
are weaker than assumed in Fig. 6 (as is likely), or if the 
characteristic ejecta velocity is a smaller fraction of the 
impact velocity (which is quite possible). Substituting the 
collision statistics developed in the run for Fig. 3 marginally 
increases the erosion times over what is shown in Fig. 6, but 
does not materially affect our conclusions. 

Therefore, unless the population of sub-km objects (which 
dominate the collision rates on comets in our model) was 
severely depleted below that predicted by the NOM and CM 
power laws, these results imply (i) that objects with radii 
— 1-2 km and smaller arc probably not mechanically primor- 
dial and (ii) that a change in the slope of the size distribution 
probably occurs for radii below —2-4 km. Depending on the 
slope structure of the primordial KD population power law, it 
may also be that the present-day disk contains far fewer com- 
ets than in the distant past. 

One factor that could stymie the collisional cascade 
among small bodies in the —30-60 AU region of Kuiper 
Disk would be a sharp cutoff in the number of small bodies. 
The recent detection of comets in the 40 AU region by 
Cochran et aL (1995) provides strong evidence that any such 
cutoff must occur below the HST detection threshold, which 
corresponds to a radius near 6 km. To test this hypothesis, 


another run was made using the fifth disk input case shown 
in Table 2. In this run, the population of bodies in the Kuiper 
Disk was fully truncated below 1 km. As shown in Fig. 7, the 
resulting survival time scales against erosion for 1 km ob- 
jects increase to much longer than the age of the solar sys- 
tem, even for eccentricities as high as 20%. As such, it can 
be concluded that the collisional cascade indicated by the 
results shown in Fig. 6 can be prevented if the Kuiper Disk 
population is somehow severely truncated below 1 km. If 
this is in fact the situation in the disk today, then it implies 
that either the number of sub-km KD bodies has always been 
severely depleted (i.e., there was a primordial size cutoff 
below 1 km), or that this condition arose through subsequent 
collisional evolution. 

Whether in fact collisions caused a depletion of sub-km 
sized objects to develop, or as may be more likely, collisions 
have created a collisional cascade to develop at sizes around 
a few km and less, two facts remain clear: First, collisional 
evolution has played a key role in shaping the population 
structure of the Kuiper Disk we observe today. And, second, 
that the signature of this collisional evolution should reveal 
itself in a distinct break in the population structure of the 
Kuiper Disk for objects with radii somewhere between ~1 
and ~6 km. 


7. SUMMARY AND DISCUSSION 

The results obtained in this paper provide strong evidence 
that collisions have been an important evolutionary mecha- 
nism in the Kuiper Disk. Indeed, because the dynamical in- 
stability time in the disk beyond —42 AU exceeds the age of 
the solar system (e.g., Levison & Duncan 1993; Duncan 
et al. 1995; Morbidelli et al. 1995), collisions appear to be 
Ae dominant evolutionary mechanism in the disk, at least 
inside 60 AU. 

The most important results obtained from the first- 
generation collision model described in this paper are as fol- 
lows. 

(1) That the total rate of collisions of smaller bodies with 
QBj-class objects is so small that there appears to be a 
dilemma in explaining how QBj’s could have grown by 
binary accretion in the disk as we know it 

(2) That present-day eccentricities in the disk preferentially 
promote erosion over accumulation for objects a few km 
in radius and smaller. 

(3) As a result, it appears that either the population of ob- 
jects smaller than ~1 km in radius was originally defi- 
cient, or the present-day population structure of the Disk 
is involved in a collisional cascade; if that later is the 
case, then many Kuiper Disk comets may not be struc- 
turally primordial. And, 

(4) That, owing to the frequency and energetics of collisions 
between several-km class and smaller bodies, a distinct 
break in the population structure of the Kuiper Disk 
likely occurs for objects with radii somewhere between 
— 1 and ~6 km. 
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FK3. 6. The time scale against erosion for 1 km radius objects at 35, 45, and 55 AU in the NOM/ CMB case, as a function of their (e), compared to the age 
of the solar system (shown as the horizontal bar). Both weaker and smaller objects erode even mm rapidly. 


Conclusion (1) is particularly important. Simply put, it 
implies that collisions appear to be too infrequent to accu- 
mulate QB r sized objects in anything approaching the age of 
the solar system. This appears to provide evidence that either 
the mass and population structure of the Kuiper Disk have 
strongly evolved over time, or that large objects like the 
QBj’s were not built via the aufbau (i.e., “building up”) 
process of binary accretion. 


Together, findings (1)— (3) strongly suggest that something 
fundamental is missing in our present state of knowledge 
about the Kuiper Disk. One possibility is that the QB| -class 
bodies were formed directly from the nebula, rather than by 
binary accretion of smaller objects. Alternatively, two possi- 
bilities based on the temporal evolution of the Disk suggest 
themselves. 

First, it may be that the mass of solids in the disk was 
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Flo. 7. Erosion time scales at 35, 45, and 55 AU in the Kuiper Disk, as a function of («), for the final collision run case shown in Table 1. In this case the 
population structure is truncated for objects with radii < 1 km. With the population truncated this way, cometary bodies suffer fewer collisions and therefore 
easily survive for longer than the age of the solar system, even for eccentricities high enough to promote erosion. 


much higher in the past than in the present A higher mass 
and therefore a higher mass density would have promoted 
faster growth of QBj bodies. The upper panel in Fig. 8 shows 
the lower-limit growth times for such a case, with 
^? disk = 123 J&q . This disk mass would be consistent with a 
continuation from 30 to 60 AU of the rather smooth surface 
mass density power law for solid material that extends from 
Jupiter to Neptune, but is today truncated at 30 AU. The 


lower panel in Fig. 8 clearly shows that “restoring the miss- 
ing mass” in the 30-60 AU zone does indeed reduce the 
lower limit to QB! growth times sufficiently. However, be- 
cause collisions between small bodies would still be ero- 
sional in a higher mass disk with such (e)' s, adding mass to 
the KD region is not (alone) sufficient to solve the QBj di- 
lemma. 

Much lower eccentricities could provide a remedy, how- 
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eccentricity of up to 3%. calculations not snown nerc aemunsumc 6 ~ * ~ LT* ^ , , ' ' . 

However, L described in Secs. 6 and 7, eccentricities below -1% arc required to permit km-class bodies to grow. The lower panel show, growth time scales 
in the same disk with a very low (*>~6X 10' 3 , which permits efficient collisional accretion from km-scale bodies upward. 


ever, by converting the collisional regime from an erosional 
state to an accretional state favoring accelerated growth. This 
is shown in the bottom panel of Fig. 8, which is a calculation 
using the same input disk as in the top panel of Fig. 8, but an 
(e) low enough to ensure efficient growth. If such low ec- 
centricities were in fact extant early in the history of the solar 


system (e.g. # before perturbations excited orbits in the disk or 
when significant nebular gas was still present), then the 
growth of larger objects would be promoted (owing in part to 
the gentler nature of collisions, and in part to the enhanced 
role of gravitational focusing at low relative velocities). 

Determining whether a higher disk mass and/or lower 
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disk eccentricities could have resulted in the growth of the 
QBj’s work requires the development of more sophisticated, 
time-dependent models that incorporate both velocity evolu- 
tion and a complete representation of the accretion process. 
We are proceeding on the latter front now. 

Before closing, however, it is useful to point out that the 
results obtained here suggest the intriguing possibility that 
the present-day Kuiper Disk shed considerable mass as it 
evolved through a more erosional stage reminiscent of the 
disks around the A stars /? Pictorus, a PsA, and a Lyr. If so. 


868 

our Kuiper Disk might be considered an older remnant of 
such a disk. 

My colleagues Don Davis, Martin Duncan, Hal Levison, 
and Glen Stewart provided useful insights during this work. 
Paolo Farineila and Eli Dwek also provided helpful com- 
ments on an early version of this manuscript This research 
was supported by the NASA Origins of Solar Systems Pro- 
gram. 
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Astronomer finds young, growing planets are easier to detect than mature ones 

San Antonio, Texas — October 23, 1996 — One of the hottest, and most 
publicly exciting areas of planetary research today is the search for planetary systems 
around other stars. Astronomers and planetary scientists want to know how 
common planetary formation is, what the range of solar system architectures is, and 
how common Earth-like bodies are in the Galaxy. 

In a significant finding. Dr. Alan Stem of the Southwest Research Institute, 
which is based in San Antonio, Texas, has found that young planets may be easier to 
detect than older ones, and that some planned groundbased telescope facilities, such 
as die Keck Interferometer in Hawaii, have the capability to observe young planets 
orbiting their parent stars by virtue of the heat the young objects give off. This 
research was recently published in The Astronomical Journal, under support from’ the 
NASA Origins of Solar System Program. 

Detecting planets around stars has long been an observational challenge 
fraught with difficulty because extrasolar planets normally only reflect light and are 
therefore intrinsically difficult to directly detect. (Typically, an Earth-like planet is a 
million or more times fainter than its parent star.) 

Dr. Stem's work examined the detectability of planets, particularly, Earth- 
sized planets, during the short but unique epoch of giant impacts that is a hallmark 
of the standard theory of planetary formation. This period is believed to have lasted 
some 30 to 60 Myr in our solar system. 

- more - 

Southwest Research Institute • 6220 Culebra Road • P.O. Drawer 28510 • San Antonio, Texas 78228-0510 • Phone (210) 684-5111 • Fax (210) 522-3547 • http://www.swn.org 



1st add — growing planets 

Sufficiently large impacts during this era, such as that believed to have 
formed the Earth's moon when a Mars-sized object struck Earth, are capable of 
turning the entire planet molten (i.e., at temperatures of 1,500 to 2,000 degrees K) 
and its atmosphere luminous, in some cases for periods of between several hundred 
to 1,000 years. While in this state, a young planet can be detected by its infrared 

radiation, which can be up to 10,000 times greater than when the same planet is cool 
later in its life. 

Stem's work found that thermally luminous Earth-sized objects can be 
detected in nearby star forming regions (which are about 125 parsecs, or almost 400 
light years away) in one to two nights of observing time. However, because even 
young planets are only sporadically heated by the truly enormous impacts needed to 
turn their surfaces molten, predictions indicate that about 250 young stars would 
have to be searched to expect to find one hot, terrestrial-sized planet. A dedicated 
observing program using, for example, 20 percent of the Keck Interferometer's 
observing time for 5 years, could find 1 to 10 such objects. These results suggest a 
new strategy for the detection of young solar systems and also offer, for the first 
time, the potential to confirm the standard model of late-stage planetary accretion, 
which involves large impacts between forming planets. 


For more information, contact: 
Dr. Alan Stem 

Geophysical, Astrophysical, and 
Planetary Sciences Department 
Southwest Research Institute 
1050 Walnut Street, Suite 429 
Boulder, CO 80302 


Tel: 303/546-9670 
Fax: 303/546-9687 

alan@everest.space.swri.edu 


### 

SwRI is an ^dependent , nonprofit , applied engineering and physical sciences 
research and development organization with a staff of 2^00 and an 
annual research volume of $243 million. 




When Oort Clouds Mingle: An Interstellar Comet Shower 


Comets roam interstellar 
space unattached to other 
stars. However, comets don't 
have to be free of their parent 
stars to pass through our 
solar system. Oort Clouds, 
the shells of comets that sur- 
round stars, can intermingle. 

Oort Clouds can be pic- 
tured as little bubbles in 
space, full of comets, with a 
diameter of roughly 20,000 
AU, about of a light-year. 
Whenever two stars with Oort 
Clouds pass within 20,000 AU 
of each other, their comet 
clouds temporarily intersect. 
Such encounters last about 


3,000 years on average. 

If half the stars in the Milky 
Way have Oort Clouds like 
our own, then about 500 
encounters with foreign Oort 
Clouds have occurred over 
the age of the solar system. 
That's about one every 10 
million years, on average. Of 
course, Oort Clouds may not 
exist around half of all stars, 
so these kinds of events 
could be much rarer. But if 
even one star in 1,000 has an 
Oort Cloud, then our solar 
system should have seen at 
least one comet shower from 
another sun. 


During each encounter, up 
to 1,000 comets from the 
intruder Oort Cloud could be 
passing within 10 AU of the 
Sun at any given time, about 
30 of which would be close 
enough to exhibit tails. Ama- 
teur astronomers would walk 
around like insomniacs from 
staying up late and looking at 
30 comets every night. The 
odds of Earth being hit by 
any of these comets, how- 
ever, are minuscule; only 
about one comet should hit 
Earth for every 300 intruder 
Oort Clouds that pass 
through our solar system. 


because the rate at which they collide with smaller 
bodies increases. Second, it also begins to dramatically 
affect the orbits of bodies coming close to, but not col- 
liding with these objects. 

It's just this action that generates the Keystone 
Cops scenes. As you can imagine, small bodies, rang- 
ing in size from rocks to Manhattans, become corks 
adrift in a rather rough sea. With many protoplanets 
growing at once, these small bodies become billiard 
balls careening among the protoplanets. 




In the shell of space between 
here and Alpha Centauri there 
may be as many as 50 trillion 
interstellar rogue comets. 

In the case of giant planets like Jupiter and Saturn, 
their gravitational effects are so strong that objects 
passing near them are frequently ejected into interstel- 
lar space. (This still occasionally happens today when 
a comet comes close to one of these behemoths). Of 
course, not all objects are ejected. Some, like Comet 
Shoemaker- Levy 9, actually do strike planets and help 
them grow. Others are scattered inward, closer to the 
Sun, where it's much harder for them to escape. Still 
others are not ejected because they don't quite get 
enough boost; these planetesimals end up in long, 
lonely orbits that can reach to distances as great ns 
tens of thousands of astronomical units (AU) from the 


Sun. This process is the one that populated the Sun's 
Oort Cloud with the icy planetesimals we call comets. 

Scott Tremaine of the Canadian Institute for Theo- 
retical Astronomy has recently shown that Jupiter and 
Saturn were so massive they tended to eject most of 
the planetesimals that came close to them. However, 
he also found that less massive Uranus and Neptune 
threw a much greater fraction of the planetesimals in 
their region of the solar system into the Oort Cloud, 
rather than onto one-way trips to galactic exile. 

Other studies, by Wing Ip of the Max-Planck-Insti- 
tute in Germany and Julio hernande/ of the University 
of Montevideo in Uruguay, indicate that Jupiter and 
Saturn probably ejected about 10 objects from the 
solar system for every one they scattered to the Oort 
Cloud. But Uranus and Neptune only ejected one 
object for every two or three they injected into the 
cloud. Overall, the process of building the giant plan- 
ets in our solar system is estimated to have injected 
several trillion (i.e., up to 3xl0 ] -) planetesimals into 
the Oort cloud, and between 10 and 100 trillion plan- 
etesimals into interstellar space. This planet-building 
stuff really is a messy business! 

Put another way, our solar system alone ejected so 
many objects into the Galaxy that the number (but cer- 
tainly not the mass) of comets in the Milky Way vastly 
exceeds the number of stars, if solar systems like our 
own are common, then this scenario has repeated 
itself billions of times, and the population of interstel- 
lar comets is impressive indeed. 

Life in the Interstellar Outback 

What is life like for interstellar comets? It isn't very 
exciting. They are subject to an eternal deep freeze 
that for all practical purposes puts them in permanent 
and very effective long-term cold storage. 

The main heat sources are the cosmic microwave 
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"background (a 3° bath from the Big Bang) and plain 
* v old starlight from the dark interstellar sky. Together, 
these two feeble radiations aren't likely to warm 
comets much above 5° or 6° C above absolute zero. At 
these cryogenic temperatures, chemical reactions are 
so slow as to be effectively nonexistent, and none of 
the common cometary surface ices that turn into gas 
when heated in the planetary region, like water, car- 
bon monoxide, nitrogen, and methane, show any 
activity. You might say deep interstellar space makes 
a very nice morgue, preserving our friendly comets 
for all eternity. 

Well, not quite. Research in the last decade or so 
has revealed that a few types of very subtle changes 
can occur on the surfaces of interstellar comets. TTte 
first breakthrough came when Bob Johnson at the Uni- 
versity of Virginia, Lou Lanzerotti of Bell Labs, and 
their co-workers showed that cosmic rays and ultravi- 
olet radiation from distant stars will penetrate the 
upper surface layers to a depth of perhaps a meter or 
two, driving out lightweight volatile molecules and 
creating micro-flaws in ice crystals. These kinds of 
radiation damage may darken and perhaps redden 
the icy cometary surfaces over billions of years. 

A few years later, work I did with Mike Shull of the 
University of Colorado showed that two other effects 
are also important. First, passing hot and massive O 
type and B type stars, and nearby supemovae explo- 
sions, will occasionally heat interstellar comets to 
comparatively balmy temperatures — perhaps 30° C 
above absolute zero (that's still about -440° F). At 
these very cold temperatures, which are ten times 
what interstellar comets normally experience, it's pos- 
sible for noble gases and a few molecules like nitrogen 
and carbon monoxide to leak out of the surface layers. 

More importantly, we also found that micro- 
impacts from smoke-sized interstellar grain particles 
will erode the surfaces of interstellar comets, perhaps 
removing their outer, radiation-damaged rind. But 
that's it. It's like the Middle Ages in Europe — only 
worse: Time goes by, but nothing cz\r changes. 


The Number of the Beasties 
The wonderful thing about interstellar comets is 
that they are a direct product of giant planet forma- 
tion. According to what we now know, if you form 
giant planets, you eject a lot of interstellar comets. The 
number of interstellar comets produced per solar sys- 
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Deep interstellar space 
makes a very nice morgue, 
preserving our friendly 
comets for all eternity. 

tern depends strongly on the spacing and masses of 
the giant planets, the mass of the disks from which the 
planets formed, the race between the rate at which 
giant planets form and begin to eject comets, and the 
length of time the disk-like planetary nebula retains its 
gas. Still, it's an exciting prospect to think that by mea- 
suring the population of nearby interstellar comets we 
can get information on the total number of solar sys- 
tems in the Milky Way's disk that have giant planets. 

But exactly how many interstellar comets should 
one expect in the Milky Way? There are two ways to 
get a rough handle on the size of this population. First, 
one could suppose all 200 billion stars in the Milky 
Way have formed Oort Clouds just like our own, each 
star ejecting some 30 trillion or so planetesimals. That 
would imply a population of roughly 6 X 10 24 inter- 
stellar comets (that's 6 trillion trillion of the beasties, 
with a total mass of about 170 million Suns). 

To get a second estimate on the number of comets, 
one can use the fact that no interstellar comet has ever 
been seen passing through the solar system. 
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This gives astronomers an upper limit on their concen- 
tration in space. From this, Zdenek Sekanina of the Jet 
Propulsion Lab calculated that on average, there is no 
more than one interstellar comet for every 1,500 cubic 
astronomical units of space near the Sun. We can use 
this number to crudely estimate the upper limit by 
multiplying the concentration of comets (Visuo, or 
0.00058 per cubic AU) and the volume of the galaxy 
(about 200 billion cubic parsecs, or 1.5xl0 27 cubic AU), 
which gives one trillion trillion (lO 24 ) comets. 

It's surprising that these two estimates are in rough 
agreement. It means that the observational data don't 
tell us very much about the population of interstellar 
comets. As much as 20 percent of the stars in the 
Galaxy could have produced Jupiter- and Saturn-mass 
planets and we still wouldn't have run across an inter- 
stellar comet in our normal comet hunting. 



It could take 1,500 years at 
the present rate of comet 
discovery to find an 
interstellar comet. 

But don't take that 20 percent number to the bank. 
Our solar system might have ejected an unusually 
large number of planetesimals, in which case the exist- 
ing constraints on the interstellar comet population 
could be consistent with jovian planets around every 
star in the whole Galaxy. Simply put, we just don't 
know enough to know — we must go and look. 

Visitors from Another Pond 

Interstellar comets would be easily distinguishable 
from solar system comets because they will pass 
through the solar system with, on average, the same 
speed that the Sun is making against the local stars. 
That's about 20 to 30 kilometers per second, compared 
to 10 km per second for comets from our own Oort 
Cloud. Because of their high velocity, interstellar 
comets follow a hyperbolic trajectory, unlike the 
parabolic trajectory coursed by solar system comets 
such as Hale-Bopp. The hyperbolic velocity of 20 to 30 
kilometers per second would make an unmistakable 
John Hancock for interstellar comet confirmation. 

Comet hunters regularly detect comets from our 
Oort Cloud, but after 250 years of comet hunting, no 
comet with a clearly interstellar orbit has been found. 
Tom McGlynn and Bob Chapman of NASA's God- 
dard Space Flight Center have estimated that it could 
take 1,500 years at the present rate of comet discovery 
to find an interstellar comet, or to prove they are rare. 
We can do far better if we try an active search. 

How would one conduct a search for interstellar 
comets? With a little ingenuity. As noted earlier, the 
best upper limit on the space density of interstellar 



comets is about one per 1,500 cubic AU. That means 
that the mean distance between interstellar comets in 
the Milky Way could be as little as 11 AU. It might 
even be a little higher near the Sun, because the Sun's 
gravity will attract them toward our direction. That's 
about the distance from the Sun to Saturn, which 
means that at any time (such as now) there should be 
one interstellar comet somewhere in the shell of space 
around the Sun as defined by Saturn's orbit. 

At a distance of 11 AU from the Sun a comet like 
Hale-Bopp would be inactive, and therefore would be 
a dark nucleus reflecting the diluted light of the dis- 
tant Sun. As seen from Earth, it would have a visual 
magnitude between 22 and 25. This is faint, but not 
too faint to discourage us; most of the recently-discov- 
ered Kuiper Disk objects have magnitudes this faint. 
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But where to search? Scouring the entire sky to 24th 
or 25th magnitude is a little much to ask. Fortunately 
we know something about where the needles are in 
this haystack. Interstellar comets will come from the 
apex of motion in the direction the Sun is moving in 
space, toward Hercules. It's a little like a meteor 
shower radiant, except for the fact that the so-called 
" meteors" are comets, and they will appear at the rate 
of perhaps one every few years. By searching a 1 
degree-wide-strip about 45 degrees in radius centered 
on the apex of solar motion about once a week, one 
can ensure that no interstellar comet will slip by. Can- 
didate objects can be followed up with subsequent 
observations after they pass through the detection 
strip in order to confirm whether or not they are in 
fact on hyperbolic trajectories from interstellar space. 


With a search strategy like this, astronomers could 
double the present-day interstellar comet detection 
limit after about 18 months. After 10 years, a ten-fold 
improvement could be achieved. By that time, it's pos- 
sible, maybe even probable, that a real, bona fide 
interstellar would have been bagged, giving us not 
only a chance to study a comet from another solar sys- 
tem, but also some solid evidence about the galaxy- 
wide frequency of solar systems with giant planets 
like our own. □ 


Planetary astronomer S. Alan Stem directs the Southwest 
Research Institute's Geophysical , Astrophysical , & Plane- 
tary Sciences research group in Boulder, Colorado. He 
wrote " Wish Upon a Star" in the December 1996 issue. 
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Chiron and the Centaurs: escapees from 
the Kuiper belt 

Alan Stern & Humberto Campins 


The Centaurs— a group of objects orbiting chaotically among the giant planets of our Solar 
System — appear to be a population transitional in size between typical short-period comets 
and the large Kuiper-belt objects that orbit beyond Neptune. They promise to reveal much about 
the origin of and interrelationships between the icy bodies of the outer Solar System. 


The outer Solar System has long appeared to be a largely empty 
place, inhabited only by the four giant planets, Pluto and a 
transient population of comets. In 1977 however, a faint and 
enigmatic object — 2060 Chiron — was discovered 1 moving on a 
moderately inclined, strongly chaotic 51-year orbit which takes it 
from just inside Saturn's orbit out almost as far as that of Uranus. 
It was not initially clear from where Chiron originated. 

Following Chiron's discovery, almost 15 years elapsed before 
other similar objects were discovered; five more have now been 
identified 1 . Based on the detection statistics implied by these 
discoveries, it has become clear that these objects belong to a 
; significant population of several hundred (or possibly several 
thousand) large icy bodies moving on relatively short-lived 
orbits between the giant planets. This new class of objects, 
known collectively as the Centaurs, are intermediate in diameter 
between typical comets (1-20 km) and small icy planets such as 
Pluto (~2300 km) and Triton (~2,700 km). Although the Cen- 
taurs are interesting in their own right, they have taken on added 
significance following the recognition that they most probably 
originated in the ancient reservoir of comets and larger objects 
located beyond the orbit of Neptune known as the Kuiper belt. 

Origin of the Centaurs 

The first clue to the origin of the Centaurs came about as a result 
of dynamical studies of Chiron's orbit. At first discovered in the 
late 1970s 3 , and forcefully reiterated in more modern calcula- 
tions 4 , Chiron’s orbit is highly unstable to perturbations by the 
giant planets. As a result, Chiron's orbital lifetime among the giant 
planets is short, leading to the conclusion that its origin was in a 
more stable reservoir, either in the asteroid belt, or beyond the 
giant planets. The discovery of a coma around Oiiron , in the 
late 1980s, indicated the presence of surface volatiles which could 
not have survived the age of the Solar System in the comparatively 
warm asteroid belt 1 ; such volatiles therefore strongly indicate that 
Chiron originated in a distant reservoir, beyond the giant planets. 

A second line of evidence relating to Chiron's origin came 
about from simulations of cometary dynamics. These studies 9 " 12 
demonstrated that the dominant dynamical class of short-period 
comets, called the Jupiter-family comets (JFC) cannot be derived 
from the classical, Oort-cloud cometary reservoir. The reason for 
this is that their characteristically low orbital inclinations cannot 
be efficiently produced by the action of planetary perturbations on 
orbits initially in the inclination-randomized (that is, nearly 
spherical) Oort cloud. 

Instead, the JFC seem to derive from a dynamically stable 
reservoir concentrated near the plane of the planetary system. 
Any such reservoir for the JFC must satisfy the criterion that the 
loss rate of objects from it be low enough that the reservoir can 
persist for the age of the Solar System. Because the dynamical 
clearing time for orbits between the planets is characteristically 
one to two orders of magnitude shorter than the age of the Solar 
System 12 " 14 , there are few regions of space that provide stable, 
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candidate source locations for the JFC. One such region is the 
zone beyond Neptune’s orbit (at 30 au, astronomical units) where 
nonlinear perturbations by Neptune and the other giant planets 
can excite orbital eccentricities on timescales comparable to the 
age of the Solar System. Once orbital eccentricities are excited 
sufficiently to cause objects to cross Neptune’s orbit, a fraction of 
these objects become temporarily trapped on Centaur-like orbits 
among the giant planets. Other such objects are dynamically 
transported by subsequent encounters with the other giant planets 
onto orbits that pass within 1-2 au of the Sun 15 , where they 
generate comae and become easily detectable. A second possible 
source region for the JFC is the slowly dynamically evaporating 
jovian Trojan clouds, whose dynamics are controlled by the 
stability of a narrow phase space surrounding the leading and 
trailing lagrangian points of Jupiter. However, recent dynamical 
simulations 16 show that the jovian Trojan clouds are not as 
effective as the so<alled trans-neptunian zone in populating the 
JFC and Centaur populations. Following historical suggestions 
dating back as far as the 1940s that a disk-shaped reservoir of 
planetesimals and other small objects might reside beyond the 
Neptune 17,11 , the trans-neptunian region has been dubbed the 
Kuiper belt, or in analogy to debris belts around the other stars, 
the Kuiper disk. 

The pivotal breakthrough concerning the reality of the hypothe- 
sized Kuiper belt came in 1992, with the discovery of a faint (R- 
band magnitude near 23), 180-km diameter object 1 * designated 
1992QB!. 1992QBj orbits the Sun in a stable, nearly circular orbit 
some 14 au beyond Neptune. In the four years since 1992QB l was 
found, over three dozen similar objects have been discovered in 
the trans-neptunian region 2 . Estimates 2 * 20 indicate that some of 
these objects have diameters approaching 400 km. Based on the 
efficiency with which such objects are being detected and their 
surface density on the sky, it has been estimated that around 
7 x 10 4 objects with diameters greater than 100 km orbit between 
30 and 50 au (ref. 2). Here we refer to these larger objects 
populating the Kuiper Belt as QB,s. 

Following the discoveries of numerous QBjS in the Kuiper belt, 
the Hubble Space Telescope (HST) was used to conduct a search 
for the much smaller, and much fainter, cometary nuclei which 
must be present in this region if it is indeed a source of the JFCs. 
Last year, Cochran et al . reported exciting evidence, near the 
limit of HSTs capabilities, for numerous objects with V-band 
magnitudes of ~28.6, corresponding to comet-like diameters of a 
few kilometres to ~10 km. This evidence corresponds to a popula- 
tion of several hundred million comet-sized objects. If this result is 
coupled with models 12 that predict the ratio of the population 
detected in the region searched by Cochran to the entire trans- 
neptunian zone, a total population is calculated of several billion 
comets in the 30-50 au region. As such, it appears that the Kuiper 
belt is indeed the source region of most JFC, as dynamical 
simulations predicted 9 u . 

Taken together, the discovery of both QB,- and comet-sized 
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TABLE 1 Orbital characteristics of the known Centaurs 


Object 

Semi- major axis 

Perihelion 

distance 

Eccentricity 

Inclination 

Present opposition 
V magnitude 

2060 Chiron 

13.70au 

8.46 au 

0.38 

25° 

15.5 

5145 Pholus 

20.30 au 

8.68 au 

0.57 

V 

17.9 

1993HA2 

24.73 au 

11.84 au 

0.52 

16° 

21.0 

1994TA 

16.82 au 

10.69 au 

0.31 

5° 

23.8 

1995DW2 

25.03 au 

18.84 au 

0.25 

4° 

21.9 

1995GO 

18.14au 

6.79 au 

0.62 

18° 

20.3 


These characteristics are taken from ref. 53. Although the basic orbital properties of the third to sixth objects 
are known, they have not yet been named because, by IAU convenbon, objects must first be observed for long 
enough to produce astrometrically reliable orbits. 


objects in the Kuiper-belt region 
indicates that the Kuiper belt 
supplies a wide size range of 
objects onto orbits in the giant- 
planet region. Based both on 
expectations resulting from the 
pianetesimal accretion codes, 
and the observational evidence 
for many more comets than 
QBjS, it appears that a power- 
law-like source population his- 
togramm exists in the Kuiper 
belt, with many more small 
bodies than QB,s. Because the 
dynamical transport process that 
brings objects from the belt to 
planet-crossing orbits is essentially independent of the mass of 
the object being transported 15 , it is expected that the population of 
objects ranging in size from Centaurs down to JFC orbiting among 
the giant planets is representative of the population of objects in 
the 30-50 au zone from which they are derived. 

Physical attributes of the Centaurs 

It is now established that the slow leakage of objects from the 
Kuiper belt due to planetary perturbations creates a population 
of objects on comparatively short-lived, planet-crossing orbits in 
the giant-planet region between 5 and 30 au from the Sun. 
Studies of the dynamical evolution of orbits dislodged from the 
Kuiper belt 15 predict a characteristic equilibrium population of 
objects on planet-crossing orbits that is ^lO" 4 of the population of 
the Ruiper-belt reservoir from which they are derived. These 
studies also predict that the median lifetime of such orbits is of the 
order of 5 x 10 7 years. Such findings imply that a population of 
~5 x 10 5 to perhaps 10 4 comets, and ~30-300 Centaur objects of 
diameter 100 km or larger, are orbiting between the giant planets. 

Chiron and its recently discovered cohort of Centaurs are thus 
now seen to be objects derived from the Kuiper belt Table 1 
summarizes the orbital attributes of the six known Centaurs; Fig. 1 
depicts the orbits of these objects and their dynamical context in 
the outer Solar System. 

As escapees from the Kuiper belt, the Centaurs are an impor- 
tant population for study. Indeed, owing to their greater proximity 
to the Sun, the brightest centaurs are some 5-7 astronomical 
magnitudes (factors of ~100 to ~600) brighter than the brightest 


Kuiper-belt objects, which enables more detailed studies of the 
Centaurs than are possible with the QB^ and Kuiper-belt comets. 
Additionally, being closer to the Sun, the Centaurs experience 
greater heating, which generates characteristic perihelion surface 
temperatures in the range ~120 to 150K (ref. 22); by contrast, 
Kuiper-belt objects probably never experience surface tempera- 
tures in excess of 60-70K. Therefore, because vapour pressure 
depends exponentially on the temperature of the ice, Centaurs are 
much more likely than Kuiper-belt objects to show sublimation- 
generated activity. Although such heating causes the surfaces of 
the Centaurs to evolve chemically and physically over long time- 
scales 23 , it also causes the surface ices to sublime, and thus reveal 
valuable insights into the nature of these objects. 

Unfortunately, although the Centaurs are brighter than Kuiper- 
belt objects, they are still faint in absolute terms, so considerable 
dedication is required to obtain physical information on them. As 
a result, comparatively little work has been done to reveal their 
compositions, colours, shapes, rotational properties and other 
attributes (Table 2). Despite the great deficits in our knowledge 
about the physical and chemical characteristics of this unique 
population, several important pieces of information are emerging. 

First, with regard to the derived sizes of the Centaurs discovered 
to date, roughly half appear to be near 60 km in diameter, but 2060 
Chiron 2124 and 5145 Pholus 23 are much larger, with ~180-km 
diameters that are comparable to typical QB,s being discovered 
in the trans-neptunian zone. Second, infrared spectroscopy and 
colour photometry have given the first clues about the surface 
compositions of these objects. The first clearly detected spectral 



FIG. 1 The orbits of the giant planets (black lines), the six known Centaurs 
(red lines) and those Kuiper-belt objects with well established orbits (green 
lines). The dot on each orbit depicts the current location of the object For 
scale, Jupiter's orbit Is approximately 10 au across. Abbreviations on the 
figure as follows: GO, 1995GO; TA, 1994TA; DW2, 1995DW2; Chiron, 
2060 Chiron; HA2, 1993HA*; Pholus, 5145 Pholus. (Figure courtesy of 
H. F. Levison) 
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TABLE 2 Physical character! sties of the known Centaurs 



Diameter* 

Geometnc 

Rotational 

Rotational 


Detected 

Object 

(km) 

albedo (in V) 

penod 

amplitude 

V - J colour} 

actrvity 

2060 Chiron 

182 ± 10 

0.11-0.20 

5.92 h 

9% 

1.13 ± 0 04 

Yes 

5145 Pholus 

185 ± 22 

0.04 ± 0.02 

9.98 h 

20% 

2.53 ± 0 06 

No 

1993HAj 

62 kmt 

0.05 



2.07 ± 0.40 

No 

1994TA 

28kmf 

0.05 




No 

1995DW2 

68 kmt 

0.05 




No 

1995GO 

60 kmt 

0.05 




No 


• The sizes of Chiron and Pholus were obtained 22 * 24 23 from thermal fluxes, with computed albedos based on their sizes and V magnitudes. Chiron ’s albedo 
Is strictly an upper limit owing to a possible, small, residual coma contribution, 
t These sizes were computed by assuming a V-band geometric albedo of 5%, a value which is commonly found for cometary surfaces. 
t Colour data are discussed in the text 27 31J2 . Not e a V-J colour of 1.116 would be identical to the Sun 47 , thereby indicating a neutrally coloured surface. 
Higher V - J colours indicate red surfaces. 


absorption feature among the Centaurs was a deep 2.25-pm 
absorption on Pholus 26 . This feature has been associated with 
light organic solids mixed with ices 273 . Importantly a 2.04-jim 
absorption features has also now been detected, both in Pholus 
and Chiron 29 , which Cruikshank et ai identify 30 in Chiron as an 
absorption band of water-ice. Third, although none of the three 
Centaurs that have been explored in the infrared have displayed 
any statistically significant colour variation with rotational phase 31 , 
they do show striking colour differences between one another 
(refs 27, 32 and Fig. 2). Indeed, whereas Chiron’s intrinsic colour 
is grey (that is, neutral) throughout the 0.3-2 J pm band, Pholus, 
which lies in a similar orbit and is similar in size, is extremely red. 
1993HAJ, though smaller and in a more distant orbit, is also very 
red compared to Chiron. How much of these differences between 
various Centaurs is due to evolutionary mechanisms (as opposed 
to intrinsic attributes) is not yet clear, but it is well established 
that long-term exposure to cosmic rays and solar ultraviolet 
radiation darkens (and initially reddens) surfaces containing 
light-weight organics, in turn creating a more complex chemical 
melange. 

Additionally, the determination of well constrained albedos for 
2060 Chiron, and more particularly 5145 Pholus (because it is not 
active), provides useful information for predicting the sizes of 
QBjS from their observed magnitudes. 

The particular value of Chiron 

Chiron is uniquely valuable among the Centaurs because of the 
long history of its sporadic outbursts. Why has such activity only 
been detected in Chiron? Possibly Chiron is dynamically younger, 
and therefore more active than the other objects. Alternatively the 
other objects may have thicker surface mantles, may be funda- 
mentally different in their composition or simply may have not 
been observed long enough to expect to detect activity. Chiron’s 



uniqueness in showing activity is perhaps the most intriguing 
observable obtained on the Centaurs so far. 

Chiron’s activity was first recognized when it suffered an out- 
burst that increased its brightness by a factor (in 1989) of just over 
two 33,34 . In 1989-91 Chiron was also observed to show a highly 
variable particulate coma and tail extending as far as 2 x 10 6 km 
(refs 35, 36), and a cloud of CN gas 37 presumably derived from the 
photodissociation of some heavier, parent molecule. When these 
various observations were made, Chiron was still more than 10 au 
from the Sun, where the solar radiation field is too weak to 
sublime water-ice, the common volatile that powers the cometary 
activity close to the Sun. Although other mechanisms remain 
plausible, the sublimation of highly volatile ices like CO, N 2 or 
CH4 (buried a short distance below the surface) were therefore 
flavoured as the source of Chiron's activity. Further evidence for 
the sublimation of such volatiles was obtained through the dis- 
covery of even more extreme activity on archival, pre -discovery 
images of Chiron obtained when it was near its aphelion at 19 au, 
and therefore far too cold to sublime anything but highly volatile 
ices like those mentioned above 3 *. The final confirmation of this 
hypothesis came in 1995, though the discovery of CO gas itself in 
Chiron’s coma 39,40 . 

The fact that Chiron’s activity was greater at its aphelion than it 
has been at any time since provides compelling evidence that its 
level of activity is not a simple function of heliocentric distance 
alone. Instead, Chiron’s activity probably involves a complex 
interaction between the level of insolation reaching its surface, 
the obliquity of its spin axis, the location of its near-surface 
volatiles and extensive surface mantling by substances (possibly 
including silicates, water-ice and carbonaceous materials) which 
do not strongly sublime that far from the Sun. 

Chiron’s strong variability and the low gas-production rate 
inferred from CN and CO gas detections in its coma provide 


FIG. 2 Colour-albedo diagram showing the visible-band albedos and 
visible-infrared (that is, V -J) colours for several Centaurs (Chiron, 
Pholus and 1993HA2), several cometary nuclei (Arend-Rigaux, Halley 
and Neujmin 1), Pluto and the one QB t (1993FW) for which applicable 
data are available. A neutrally coloured surface would have the V - J colour 
of the Sun, 1.116 (ref. 54). Note that the Centaur Pholus is quite red; In 
fact, it Is far redder than any other object in the Solar System. The V - J 
colour for 1993FW is an upper limit V - J for Pluto was obtained form D. P. 
Cruikshank (personal communication; the V-albedo error bars for Pluto 
represent its intrinsic rotational lightcurve variation. 
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strong evidence that Chiron’s activity is generated by localized 
sources covering <1% of the surface. The case supporting highly 
localized vents or jets on the surface is further supported by short- 
term brightness fluctuations in Chiron’s coma*' * 1 11 that occur on 
timescales consistent with clouds of material being ejected onto 
suborbital trajectories, and by the detection of complex opacity 
structures in Chiron’s coma during two recently observed stellar 
occultations by Chiron 24,43 . It has been pointed out 44 that these 
vents or jets may resemble the geysers detected on the surface of 
Neptune’s large, captured satellite, Triton. 

Chiron’s low gravity creates a situation in which its escape speed 
(~10 2 ms -1 ) is comparable to both the thermal velocity of sub- 
liming gas (~2x l(r ms" 1 ), and the estimated muzzle velocity of 
Triton-like geysers 45 (of the order of 40-300 ms -1 ). As a result, 
some of the gas and entrained particulates ejected from Chiron’s 
surface would be deposited into high suborbital trajectories; much 
would also escape. Modellers have only begun to explore the 
range of interesting physical phenomena likely to result in this 
intermediate regime between freely escaping cometary comae and 
strongly bound planetary atmospheres 46 . Among these is the 
distinct possibility that Chiron’s neutral colour and comparatively 
high albedo are the direct result of its activity, which probably 
causes a thin veneer of icy particulates to rain back onto the 
surface from suborbital trajectories. 

An emerging view 

We are witnessing a revolution in our understanding of the 
content and architecture of the outer Solar System. Whereas a 
decade ago, the outer Solar System seemed to consist only of the 
outer planets, the Oort-cloud comets and the then-rogue object 
Chiron, we now see revealed both the teeming Kuiper belt and its 
progeny, the comets and Centaur-sized objects orbiting among the 
outer planets. As a result, we have come to recognize that the 
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outer Solar System is littered with icy objects intermediate in size 
between comets and the giant planets. 

We have also learned that the early stages of planetary forma- 
tion, with widespread growth from planetesimals to objects with 
diameters of several hundreds of kilometres, provide concrete 
evidence for an ancient era of planet-building in the Kuiper-belt 
region 47,48 . For some reason (probably involving the role of 
Neptune that excited orbital eccentricities that were not condu- 
cive to further growth), the era of accretion in the Kuiper belt was 
prematurely trunctated at a stage where intermediate-sized 
objects had formed 49 . The strong circumstantial evidence for the 
early formation of numerous objects in the 1,000-km class, of 
which Pluto and Triton are apparently the sole extant remnants 
within observational reach 30 ’ 51 , further supports the case for 
initially strong but eventually arrested planetary accretion in the 
Kuiper-belt region 4 * 42 . As such, the Kuiper belt has become one of 
the most important regions in the Solar System for studies of 
planetary origins. The Centaurs and QB,s therefore represent a 
valuable, relic population of icy objects whose growth was arrested 
at a fascinating, intermediate stage between comets and small 
planets. 

The Centaurs also serve as bright proxies for distant comets, as 
laboratories for studying surface processes occurring on comets, 
Triton and perhaps Pluto, and as nearby proxies for the QBjS and 
other intermediate-scale bodies that bridge the size gap between 
comets, Pluto and Triton. As such, they hold special promise for 
understanding the origin and interrelationships among the icy 
bodies of the outer Solar System. □ 
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